Stability of low-friction surface sliding of nanocrystals with rectangular symmetry and 

application to W on NaF(OOl) 
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We investigate the stability of low-friction sliding of nanocrystal with rectangular atomic arrange- 
ment on rectangular lattices, for which analytical results can be obtained. We find that several 
incommensurate periodic orbits exist and are stable against thermal fluctuations and other pertur- 
bations. As incommensurate orientations lead to low corrugation, and therefore low friction, such 
incommensurate periodic orbits are interesting for the study of nanotribology. The analytical results 
compare very well with simulations of W nanocrystals on NaF(OOl). The geometry and high typical 
corrugation of substrates with square lattices increase the robustness compared to typical hexagonal 
lattices, such as graphite. 
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I. INTRODUCTION 

Friction at the atomic scale is currently actively stud- 
ied [l| . One of the goals of this research is to understand 
whether extremely low friction can be obtained by an 
appropriate choice of the sliding conditions. Commensu- 
rability between the sliding lattices is one of the elements 
that determine friction. For a purely incommensurate in- 
finite contact, theoretical arguments suggest that static 
friction should vanish 2]. However, very low friction 
has been measured also for finite incommensurate con- 
tacts at very low velocities [3j, |4j and this effect has been 
called superlubricity 0-0] • The atomic force microscope 
(AFM) study 0, Ej found that the sliding of graphite 
flakes on graphite can occur with very low friction, de- 
pending strongly on the relative orientation. At the same 
time, rotation of the flake can lead to a rapid increase of 
friction and stick-slip motion, corresponding to a locking 
into a commensurate orientation. States of very different 
friction have also been observed to coexist in the sliding 
of nanoparticles [8| and have been attributed to contam- 
ination or amorphous surfaces. 

This paper examines theoretically the sliding of 
nanocrystals and substrates which both have different, 
but regular rectangular lattices. By means of a simple, 
analytically soluble model, we show that stable orien- 
tations exist for any size nanocrystal, and derive some 
general properties, showing how the stability depends on 
the scan line. The stable orientations are independent of 
the velocity and nearly independent of the corrugation, 
but depend only on the geometry of the substrate and 
contact layer of the nanocrystal. We can also estimate 
the energy barriers necessary to rotate the nanocrystal 
from a given orientation to another. 

We apply the results to W nanocrystals on an 
NaF(OOl) substrate. NaF has been studied extensively 
as a substrate in the context of nanotribology 
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and W is commonly used for AFM tips (see, for instance 
Ref. 0). Both materials have a bec lattice structure and 
the (001) surfaces in contact can thus be described by one 
finite and one infinite square lattice. The lattice parame- 
ters of the two materials are very similar (0.31585 nm and 
0.32668 nm respectively), which ensures the existence 
of nearly commensurate orientations for small enough 
nanocrystal. For W, the corrugation on NaF (001) is of 
the order of 1 eV, a value ty pica l for many substrate 
systems (see, for instance Ref. [lOj)- 

In this paper, we apply a model previously introduced 
to study the dynamics of hexagonal flakes on hexago- 
nal lattices @ to rectangular contact layers on substrates 
with rectangular lattices. One of the important results 
is that the inherent lack of robustness of the low friction 
motion found in the hexagonal graphite system for par- 
ticular scan lines is not present for rectangular lattices, 
which opens the possibility for practical applications of 
low- friction sliding. The analytical results of the simpli- 
fied model are supported by our numerical simulations of 
W nanocrystal on NaF. 

In Sec. |TT] we introduce the notation while briefly re- 
viewing the model used in Ref. [6] to describe the relation 
between the rotational dynamics and friction. In Sec. lIIIl 
this model is applied to a general system consisting of a 
rectangular lattice and a nanocrystal with a rectangu- 
lar contact layer. The analytical results are compared 
to numerical simulations of W nanocrystal on NaF(OOl) 
in Sec. IIVI and finally the implications of our results are 
discussed in Sec. [Vj 



II. THE SIMPLIFIED MODEL OF FRICTION 
AND ROTATION 

We wish to study the dynamics of small rigid nanocrys- 
tals of atoms arranged in a regular lattice on a regular 
substrate, as sketched in Fig. [TJ The contact layer con- 
sists of one atomic plane and lies in the x—y plane parallel 
to the periodic substrate. For matched lattice parame- 
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ters, by changing the orientation 4> of the nanocrystal on 
the substrate, the contact layer can be either commensu- 
rate or incommensurate with the substrate atoms. The 
incommensurate states which produce low friction can be 
destroyed by rotations around the z axis that lead to a 
locking in a commensurate orientation 

Each atom of the contact layer is subjected to a pe- 
riodic potential due to the substrate and to an external 
load force which may be applied to the nanocrystal. Ad- 
ditionally, the contact layer atoms experience a friction 
force due to dissipation of kinetic energy into phonon 
modes of the substrate crystal. The remaining nanocrys- 
tal atoms are further away from the substrate and we 
may therefore assume that they couple only to the other 
nanocrystal atoms. The AFM cantilever or a similar de- 
vice used to investigate friction can be modelled by cou- 
pling the centre of mass of the nanocrystal harmonically 
with a spring with force constant c to a support moving 
at constant velocity v s . The friction is given by the av- 
erage value of the lateral force F s , exerted by the spring 
on the nanocrystal. 

In three dimensions, a rigid nanocrystal is thus left 
with only 6 degrees of freedom: the coordinates of the 
centre of mass, and the orientation. Despite this, the sys- 
tem is still too complicated to perform the stability anal- 
ysis analytically. In Ref. [6| a simple model was proposed 
to describe the rotation of hexagonal graphite flakes on 
hexagonal graphite lattices and its consequences for low- 
friction sliding. We briefly summarise this model and 
some of the results of Ref. [6| here in a general context. 

A possible simplification of the system is suggested by 
the nature of the dynamics. We are interested in the ro- 
tation around the z axis, which affects the commensura- 
bility, and therefore the friction. Hence, the two relevant 
degrees of freedom are the position of the centre of mass 
along the scan line, x, and the orientation, <fr. Conse- 
quently, instead of the full substrate potential V(x, y, z), 
we introduce a potential V(x,(j>), one example of which 
is given in Fig. [5J where one can see that the corrugation 
felt by the nanocrystal as a whole decays quickly away 
from the commensurate orientations cf> — 0,90°. Such a 
model is fully described by the initial support position 
x®, support velocity v 3 , mass M and moment of inertia 
I of the nanocrystal, the effective potential V(x, 4>), and 



FIG. 2. An example of the shape of V(x, (f>) as obtained 
from the potential V\(x,y) in Eq. (|15|) at constant y and z, 
summed over the atoms of a nanocrystal with square con- 
tact layer of 4 x 4 unit cells (and therefore 5x5 atoms) with 
lattice parameter 0.31585 nm on a square lattice with corru- 
gation Vx = 1 eV and V2 = V3 = Vi/2 and lattice parameter 
0.32668 nm. These parameters correspond to W(001) and 
NaF(OOl). 



the viscous friction coefficient 7 of a contact layer atom 
on the substrate. Altogether, this model is similar to 
the Tomlinson model[9j, but with one additional degree 
of freedom which accounts for the rotational dynamics. 
The details of the substrate and contact layer lattice ge- 
ometry, choice of scan line, and the applied external load 
force are accounted for in V(x, cj>) [6|. 

We write the equations of motion of the simplified sys- 
tem as a dynamical system of first-order differential equa- 
tions for the position x, velocity v x , orientation <f>, and 
angular velocity uj, 



Mv 



x = v x , 

dV{x,ct)) 



dx 



(1) 

c(x — tv s — x s o) — 'yMvx , (2) 



4> = u , 

loj = — ■ylu! 



(3) 
(4) 



Some general properties of the potential energy V(x, <fi) 
can be derived from the symmetries of the substrate lat- 
tice and contact layer. The translation symmetry of the 
substrate dictates that V(x, </>) must be periodic in x. A 
good representation of V(x, 4>) is therefore given by 



V(x, <p) = U(4>) + W(4>) cos 



/2nx 
W 



(5) 



where U(<j>) and W(4>) are both smooth functions that 
represent the average value of the potential energy and 
the amplitude of the modulation respectively. 

Furthermore, rotational symmetries of the substrate 
and contact layer lead to rotational symmetries of V(x, (f) 
of the form 



V{x,(f>) = V{ X ,(j>; 



symmetry 



(6) 
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where ^symmetry is an angle of rotation under which the 
contact layer is symmetric. 

These symmetries of the potential V(x, (f>) further im- 
ply that 



U(tf>) 

w(<p)-- 



- U (Asymmetry 
(^symmetry 



(7) 
(8) 



For the case we consider here, square lattices, we show 
later in Sec. MI that also 



U(<p) = U(- 
W((f>) = W(- 



(9) 
(10) 



This also holds for some combinations of parameters of 
rectangular lattices. In turn, these equations imply that 
U and W have extrema in </> = <fio =0, </> S ymmetry/2. Since 
the torque, given by Eq. (Q}, vanishes for oj = 0, 4> = <f>o 
these conditions define a two-dimensional invariant man- 
ifold of the dynamics. If the nanocrystal is on the man- 
ifold, it will remain there, and keep its orientation. In 
general, there may be more invariant manifolds at other 
orientations, if U and W have additional extrema which 
coincide. 

If the orientation of an invariant manifold is incom- 
mensurate, the corrugation is relatively small (see Fig. ^ 
and the friction force on the nanocrystal is small as well. 
Such incommensurate invariant manifolds are therefore 
interesting for the study of nanotribology. 

In Fig. [3J the consequences for friction are shown for a 
typical system consisting of a W nanocrystal with (001) 
face on NaF(OOl) as described by simulations of the full 
two-dimensional system with the single-atom interaction 
potential of Eq. (fT5j) . The orientation of the nanocrystal 
remains nearly constant and the motion in the x direc- 
tion exhibits stick-slip behaviour which is relatively mild 
for incommensurate orientations and corresponds to a 
relatively low average friction force. 

The existence of incommensurate invariant manifolds 
does not necessarily mean that low friction can be ob- 
served under experimental conditions. The invariant 
manifolds may not be stable, namely, the nanocrystal 
may rotate away from it if it is at an orientation which 
deviates slightly from the invariant manifold. We must 
therefore investigate the stability and robustness of the 
invariant manifolds. 

Consider a general potential V(x, </>) which has an in- 
variant manifold at cb = d>a, i.e. 



dV{x,4>) 







(11) 



for all x. The growth rates of perturbations in <j) and lo 
are equal to the associated Lyapunov exponents, which 
are 



A± 



2 ' 2\ 



d 2 V(x,<p) 



(12) 
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FIG. 3. Typical simulated trajectories of a W nanocrystal 
containing 4x4x4 unit cells (contact layer 5x5 atoms) on 
NaF(OOl) in the absence of thermal fluctuations, at y — a/2, 
being pulled in the x direction. The orientation remains 
nearly constant (a) and shows only small periodic fluctuations 
with period a/v s . Both the commensurate and incommensu- 
rate orientations display stick-slip behaviour in the position 
(b), but the friction (c) is significantly higher in the commen- 
surate state. Further details of the simulations are described 
in Sec. IW] 



where () t denotes the time average of a quantity on a typ- 
ical trajectory on the invariant manifold. If these growth 
rates are smaller than zero, the invariant manifold is sta- 
ble. This is the case if the time average of the potential 
energy is at a local minimum, i. e., 



d 2 V(x, 



) >°- 

4>=4>o I t 

Using Eq. §5§ , Eq. (fT3"]) can be rewritten to read 



(13) 



d 2 U((f>) d 2 W(<f>) 



(¥) 



> o . 

(14) 



0=00 , 



The stability is thus determined by the functions U and 
W, and how much time the particle spends near the min- 
ima of the potential, where the cosine is negative. 

In stick-slip motion, the particle spends most of its 
time in the minima of the potential, i.e. where the cosine 
is smaller than zero (see Fig- EJ) ■ In extreme cases, (cos)t 
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FIG. 4. A top view of a general rectangular lattice (open 
circles) and contact layer (filled circles) with mismatch angle 
4>, lattice parameters ai , a 2 , 61 , hi , contact layer size m,n, and 
the potential energy of a contact layer atom on the substrate. 
If a contact layer atom lies on top of a substrate atom, its 
potential energy is Vi. If it lies directly between an atom and 
its nearest neighbour in the x or y direction, it has potential 
energy V% or V3 respectively. If it lies in the centre of a rect- 
angle, at equal distance from four substrate atoms, without 
loss of generality, we may set the potential energy to 0. The 
origin of the coordinate system is chosen to lie on top of a 
substrate atom. 



may be almost equal to — 1 . If the motion is truly super- 
lubric, then the particle spends about the same time in 
the minima as it does in the maxima. If the motion is 
nearly superlubric, then the particle spends most of its 
time in the minima. Hence, for realistic cases, (cos)t < 0. 

If the offset of the potential, C/(0), has a minimum at 
(j)Q it contributes positively towards the stability. Simi- 
larly, if the amplitude W((j>) is at a maximum at <fio the 
stability is enhanced, because the second derivative is 
multiplied (cos) t, which is a negative number. A mini- 
mum of U and maximum of W therefore always lead to 
stability, whereas a maximum of U and minimum of W 
always leads to instability. If both are at a maximum, or 
both are at a minimum at c/)q, then the stability is not 
directly obvious. 

The representation of the potential used in this sec- 
tion not only gives qualitative understanding, but also 
allows for quantitative predictions about the dynamics, 
as is shown by comparison to full numerical simulations 
in Sec. [iVl 



III. U AND W FOR RECTANGULAR LATTICES 



The functions U(<p) and W((f>) determine the existence 
and stability of the low-friction states. We will now inves- 
tigate these functions for a range of commonly occurring 
substrate and lattice combinations with rectangular unit 
cells. 

Consider a two-dimensional substrate consisting of a 
rectangular lattice, see Fig. @] The potential energy of 
an atom of the contact layer at position (X, Y) due to the 
presence of the substrate can be written in the general 



form 

V A (X,Y) = 

V1-V2 
2 

, V!-V 2 



V2 + V3 
2 

ai 



V 3 



V!-V 3 



a-2 



2tt- 



X 
a 1 



2tt- 



Y 

(12 



1L5.) 



Because the potential is periodic in X with period 01, 
this can be rewritten in a form similar to Eq. ([5]), 



V A (X,Y) = U A {Y) + Rc 



W A {Y)exp 2iri 



X 
a 1 



(16) 



with Ua(Y) and Wa(Y) the average potential energy and 
corrugation of a single atom travelling in the x direction 
at constant y — Y. They can be determined from the 
potential Va(x,v) in Eq. ([T51) . 



Ua(Y) = ~V a (0,Y) + ±Va(y,Y) 

Vi + v 2 + v 3 



4 



Y 



2tt- 

0,2 

Wa{Y) = \va{Q,Y)- 1 -Va{^,Y) 



V1-V2 + v 3 
4 

, vi-y 2 



v 3 



2tt- 



Y 



(l 2 



(17) 

(18) 
(19) 

(20) 



The nanocrystal is rigid, rectangular and the contact 
layer consists ofmxn atoms. The total potential energy 
of the contact layer depends on the position (x, y) of the 
centre of mass and orientation (f>, and can be expressed 
as a sum over the potential energies of the individual 
contact atoms, 



V(x,y,(f>) 



j2J2 v ^ x ^Y ]k ) 

j=i k=i 



(21) 



where (Xj^, Y^) is the position of the atom in the j-th 
column and k-th row of the contact layer. For a rectan- 
gular three-dimensional nanocrystal the centre of mass is 
directly above the centre of the contact layer, i. e., 



Y jk 



COS ( 

sin c 



x 

y 

— sku 
cos 6 



bi (j - 
62 (k 



m+l 
2 

n+1 



(22) 



Let us write V(x, y) in a form similar to that of Eq. ([5]), 



V(x,y,</>) = U y {<i))+Rc 



.(23) 



where U y and W y are the average potential energy and 
corrugation of a contact layer at orientation <fi moving on 
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a scan line along the x direction for a specific y. They can be obtained from Ua and Wa as, 

m n 

3 = 1 k=l 

m n j 

j=i fc=i ^ 



ai 



(24) 
.(25) 



Substituting Eq. (|2"2l and explicitly performing the sums 
yields 



V 1 + V 2 + V 3 V!+V 2 -V 3 / y ^ sm^u^cos^sm^TO-sm 
Uylo) = mn 1 cos 2tt — 



sin cos 0^ sin (^rm^ sin </>^ 



sin 1 7T^- cos J sin I ir^ sm < 



Wy(4>) = 



Vx-V 2 + v 3 



Sin ( 7TTO — 



cos m sm 



sin ( 7m— sir 



sin ^tt^- cos <fij sin ^tt^- sin (j> 



V1-V2-V3 / y 

■ COS I ZTT 

a-2 



sin 


vnb (s*± - SB-tX 


sin 


\ ai a 2 J 


sin 


\ a>2 ffii 1 


sin 


1 ^ Ol £12/ 





(26) 



(27) 



where we have used that Ylf=o exp{2ix[Z — (d — 1) /2]} = 
sin(da;)/ sin(x) 

We consider the shape of these two functions U and 
W for the important case of square lattices in the next 
section. 



A. Square contacts on square lattices 

Let us consider square contact layers (to = n = d) of 
atoms arranged in a square lattice (61 = 62 = b) on a 



r 



substrate with a square lattice (ai = a% = a, V2 = V3). 
For this case Eqs. (|26l) and (|27|) become 



with 



U y (<l>)=d 2 



cos 



(2tt^) u#) , 



W y (4>) = iu((/>) + cos (2tt^) wi(0) , 



(28) 
(29) 



w(4>) 



V\ sin (nd^ cos 0) sin (7r<i^ sin c/>) 
4 sin (71"- cos 4>) sin (tt^ sin 0) 

Vi — 2V2 sin [nd^ (cos ^ — sin 0)] sin [7rd^ (cos + sin < 
4 sin [7r| (cos — sin </>)] sin [71-^ (cos </> + sin (j)' 

I 



(30) 
(31) 



These equations are symmetrical under the transforma- 
tion cj) — » — <\>. The second term on the right hand side of 
Eq. P5]). iui(0), is proportional to (Vi - 2V 2 )/^ which 
is generally small, and in many cases taken to be 
for instance, the substrate potentials used in Refs. 
We therefore initially neglect the second term and focus 
on the first term and write 



The function w(<j)) then completely determines U and W 
for all scan lines, and therefore the invariant manifolds 
as well as their stability. As the first term of U y ((j>) is 
constant, the extrema of U y coincide approximately with 
the maxima of | W y | , leading to invariant manifolds at the 
corresponding orientations, 



W y ^)=w{cj>) 



(32) 



dw{4>) 



= . 



(33) 
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FIG. 5. The function w(4>) which determines U y and W y 
(see Eqs. (|28[) and (|32[) . and thus the invariant manifolds 
and their stability, is plotted for various values of d and 
b/a = 0.96685, corresponding to a W nanocrystal with its 
(001) face on NaF(OOl). The number of extrema increases 
with d. The maximum size of the extrema for d — > oo is also 
plotted. In general, if d is even, w(<f)) has a minimum at 0, 
and if d is odd a maximum. For d = 5, besides the maxima at 
4> = 0° and (f> = 45°, w(tj)) has minima at 4> ~ 16.95° , 36.25° , 
and a maximum at ~ 28.32°. For d — 10, besides the 
minimum at <j> = 0° and the maximum = 45°, w(cj>) has 
minima at cj> « 14.33° , 22.98° , 32.56° , 40.86° and maxima at 
~ 8.39°, 19.51°, 27.69°, 36.34°. Each of the extrema corre- 
sponds to an invariant manifold. 



In Fig. [5] examples of w{<j>) are plotted for W on NaF(OOl) 
for several values of d. 

The extrema of U and W give the invariant manifolds. 
One of the invariant manifolds is at <j> — 0, while the 
others are at incommensurate orientations, which do not 
depend on the scan line y. Additionally, the nodes in 
the amplitude {W y (4>) = 0) all occur at the same value 
of U y {4>) = d 2 {Vi + 2V 2 )/4. The number of such nodes 
can be determined from the number of zeros of the nu- 
merator, and in the range of G (0°, 45°) for a = b, this 
number is in general equal to d — 1, except when the ze- 
ros coincide, which occurs for instance when d 2 can be 
written as the sum of the squares of two integers. 

As w(4>) is independent of the position y of the scan 
line, the orientations at which invariant manifolds oc- 
cur are independent of the scan line while their stability 
is controlled entirely by the prefactor cos(27ry/a) in U y . 
The requirement of stability, Eq. (fT4|) . implies that each 
minimum of U y ((j>) leads to a stable invariant manifold, 
which in turn means that there are at least (d — l)/2 
stable incommensurate periodic orbits if a = b. If the 
scan line changes and the prefactor cos(2ny / a) becomes 
small, the amplitude term of W in Eq. ([T5]) becomes 
dominant, and the invariant manifolds at the maxima 
of U become stable as well. If the prefactor cos(27rj//a) 
changes sign, U y changes sign, so that the maxima of U 
becoming minima and the minima maxima. The invari- 
ant manifolds which are stable for cos(27ry/a) > then 
become unstable for cos(27ry/a) < 0, and those which 



were unstable become stable. Depending on the sign and 
size of cos(27ry/a), there are therefore generally between 
(d— 1)/2 and d— 1 stable incommensurate periodic orbits. 
In cases of very strong stick-slip behaviour, the centre of 
mass spends most of the time in the potential minimum 
and (cos)t is almost equal to —1. In this case, the invari- 
ant manifold is stable for nearly all scan lines. 

When d becomes large, the numerator in w((j)) oscil- 
lates rapidly. However, the denominator is independent 
of d and determines the size of the maxima and minima 
(see Pig. [5]), 



<(<f>) = 



1 



V a 



(34) 



This function therefore estimates the robustness of the 
incommensurate orientations, i. e. the typical energy bar- 
rier that must be overcome to rotate from one to another 
nearby stable orientation, 



AE(<j>) 



w n 



(35) 



In the case of the graphite flakes on graphite, it was 
found that the incommensurate orientations are not very 
robust 6]. However, graphite has a very small corru- 
gation (around 25 meV, comparable to k^T at room 
temperature) compared to NaF(OOT), which has a cor- 
rugation of around 1 eV. The typical energy barrier for 
W on NaF(OOl) is around 0.4 eV even at 45°. Conse- 
quently, the low-friction incommensurate states of this 
system are robust against thermal fluctuations at room 
temperature. Additionally, due to the much larger mass 
of the W atoms and large M and /, the incommensurate 
states are also more robust against low support velocities 
than in graphite. 



B. A note on hexagonal lattices 

The results presented here can be applied to determine 
U and W for a graphite flake subject to the potential 
of Ref. [H, which consists of a sum of two rectangular 
lattices, while a hexagonal flake can be written as the 
sum of several rectangular flakes. 

Similar behaviour, including the existence of the tran- 
sition scan line where cos(27ry/a) = was observed 
for hexagonal flakes of different sizes on a hexa gona l 
lattice [q. The diameter for hexagonal flakes is 2WN/6, 
and the total number of invariant manifolds is 2\J N/6+1. 
It should be noted, however, that in hexagonal lattices, 
due to the symmetry, the corrugation vanishes for some 
scan lines, whereas this does not generally happen for 
square lattices. This low W y in hexagonal lattices tends 
to occur at scan lines where the minimum in U y at 
the incommensurate orientation is not very pronounced, 
causing the incommensurate states to become only very 
weakly stable @ and not robust. 
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IV. W NANOCRYSTALS ON NAF(OOl) 

The results of the previous section are applicable to 
many systems that present square lattices, like the (001) 
surfaces of fee and bec materials. In this work, as a rep- 
resentative system, we have chosen a substrate which has 
been widely studied experimentally, NaF(OOl) llj. Due 
to the bec structure of the NaF, the surface atoms of the 
substrate are arranged in a square lattice with lattice 
parameter 0.32668 nm. Additionally, the corrugation of 
NaF, around 1 eV, is typical for many substrate systems. 
Unless otherwise mentioned, in this work corrugation pa- 
rameters of V\ = 2V2 = 2V3 = 1 eV will be used. 

For the nanocrystal, we have used W, which is also 
common in experiments as a tip material. W has a bec 
structure and the (001) contact layer consists of a square 
lattice, with lattice parameter 0.31585 nm. As the lattice 
parameters of the contact layer and substrate are very 
similar, there exist orientations at which the two lattices 
are nearly commensurate, as long as the contact layer is 
small. 

This system is thus extremely suitable for studying 
the effects of commensurability on the friction of small 
nanocrystals. As the W atoms are quite heavy, the to- 
tal mass and moment of inertia are large, and the mo- 
tion of the nanocrystal should behave approximately one- 
dimensionally and the simple model described in Sees. HTl 
and lllll is applicable. The function w(4>) (Eg. [31"]). plotted 
for several cases in Fig. [3] determines the stable orienta- 
tions. These orientations are mentioned in the caption 
for d — 5 and d = 10. 



A. simulations 

We have performed molecular dynamics simulations of 
rigid W nanocrystals on NaF(OOl) based on the full po- 
tential Va(X, Y) in Eq. (|15p. The nanocrystals were cou- 
pled harmonically (spring constant c = 10 N/m) to a sup- 
port moving at constant velocity v s = 32.668 mm/s (or 
0.1 a/ps) in the positive x direction at constant y = y s . 
Note that this is several orders of magnitude larger than 
the velocities typically used in AFM experiments, which 
are of the order of 1 fi/s. These values were chosen for 
computational reasons, but the dynamics of the system 
at lower velocities are represented well, as the time scales 
associated with the interaction between the substrate and 
nanocrystal are still much shorter than the time it takes 
the support to traverse one unit cell. For the viscous 
friction parameter of the substrate we have chosen the 
typical value of 7 = 1/ps. For simplicity we have re- 
stricted ourselves to cubic nanocrystals. 

In Fig. [21 several examples are shown of simulated tra- 
jectories for d = 5. The orientations are nearly constant, 
but may fluctuate periodically with a period equal to the 
time it takes the support to move one lattice spacing. 
These fluctuations are due to the motion in the y direc- 
tion. 
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FIG. 6. A bifurcation diagram for a cubic W nanocrystal 
with a 5 x 5 contact layer on NaF(OOl) at support velocity 
v s = 32.668 mm/s. Simulations were run starting from initial 
conditions x — —a/4,y = y B ,v x = v B ,v y = 0,u = and 
a range of initial orientations between and 45°, with 0.5° 
intervals. The final angle is plotted in black. The set of initial 
angles which converge to the various periodic orbits, i.e. cross 
sections of the basins of attraction, are plotted in different 
colours for each stable orientation. Apart from <j> = 0°,45°, 
stable orientations occur at cj> ~ 16.95°, 28.32°, 36.25°, which 
correspond exactly to the extrema of w(<f>) for this system, 
shown in Fig. [S] 
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FIG. 7. A bifurcation diagram and basins of attrac- 
tion of a system similar to the one of Fig. [6] apart 
from the size of the W nanocrystal, which is still cu- 
bic, but has a 10 x 10 atom contact layer. Apart 
from <j> — 0° , 45° , stable orientations occur at <f> « 
8.39° , 14.33° , 19.51° , 22.98° , 27.69° , 32.56° , 36.34° , 40.86° . 



In Fig. [5] the results of 11739 simulations of the same 
system for different initial conditions and y s are shown. 
The final angles is plotted as a function of the scan line 
and initial angle, demonstrating the existence and stabil- 
ity of the periodic orbits at specific orientations. These 
orientations correspond exactly to the extrema of w{4>) 
for this system, shown in Fig. [5] 

The sets of initial angles which converge towards these 
stable orientations are also indicated in Fig. HI It can be 
seen that the orientations at 16.95° and 36.35° are more 
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FIG. 8. The functions w(4>) and wi((f>) are plotted for b/a — 
0.96685, corresponding to a system of W(001) on NaF(OOl) 
for d — 5 and V2 = Vi/3. w(<f>) is the same as in Fig. 
but W y ((j)) is no longer equal to it. Besides the maximum 
at (j> = 0° and minimum at cfi = 45°, W y (4>) has minima at 
<l) sa 11.73°, 23.83°, and maxima at <p « 17.63°, 32.39°. 

stable at y s = 0, as expected from minima of w((f>), while 
0°, 28.32°, and 45°, which correspond to maxima, are, as 
expected, more stable at y s = a/2. 

In Fig. [7J the same plot is repeated for d = 10. There 
are more stable orientations, and these correspond ex- 
actly to the extrema of w(<fi) for d = 10, shown in Fig. [5] 
Their behaviour near y s — and y s = a/2 is also exactly 
as expected from w(4>). 

The orbit at (f> = 36.35° in Fig. |S] breaks down for 
some < y s < a/4. For these scan lines the theoretical 
approach described in this paper predicts that this ori- 
entation is stable. However, it is impossible to explore 
all possible sets of initial conditions in simulations, and 
thus the bifurcation diagrams are necessarily incomplete. 
At these scan lines, the range of initial conditions chosen 
for the simulations does not intersect with the basin of 
attraction of the stable periodic orbit at <fi = 36.35°. 

The comparison of these numerical simulations to the 
predictions of the analytical model is exceedingly good, 
including the important feature that several incommen- 
surate orbits are stable for all scan lines. 



B. The case V 2 = V 3 # Vi/2 

Though in most cases of square lattices V2 = Vi/2 
(see, for instance Ref. |l0|), this may not always be true. 
We therefore briefly consider the implications of V2 — 
V3 7^ Vi/2, where, in Eq. (|29p . the second term can no 
longer be neglected. In Fig. [8] w{<j)) and w\ (</>), which 
determine U y {4>) and W y {4>) through Eqs. (J^J) and ((21))) 
are plotted. These two functions together determine the 
invariant manifolds and their stability. The extrema of U 
and W no longer coincide exactly, and so the location of 
the invariant manifolds is much less trivial than before. 

In Fig. 13 the plot of Fig. [5] is repeated for V2 = Vi/3. 
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FIG. 9. The plot of Fig. |6] repeated with Vi = Vi/3. 
Besides at <f> = 0° , 45° , stable orientations occur at <f> ~ 
16.89°, 27.9°, 36.59°. 



As both w(4>) and W y (4>) still have an extremum near 
4> « 17°, the orbit at that orientation persists, but is 
changed slightly. The orbits near 28° and 36° also still 
exist, in distorted form, because the second term in W y , 
determined by wi((f>) cx (V\ — V2VV2, is relatively small 
compared to U y {4>) 1 so that U y ((j>) is the dominant contri- 
bution to the stability, regardless of the scan line. This is 
in general the case, unless V2 deviates very strongly from 
Vi/2. 



V. CONCLUSIONS 

We have analytically investigated the driven nonlinear 
dynamics of general rectangular nanocrystals on rectan- 
gular lattices and their relation to sliding friction. We 
have formulated an approximate analytical model, which 
gives the conditions for the existence and stability of in- 
commensurate sliding states with low friction. We show 
that the number of incommensurate orbits grows linearly 
with the diameter of the nanocrystal and that for realis- 
tic systems, several states can be robust against thermal 
fluctuations, change of scan line, and driving velocity. 
The geometry and high typical corrugation of substrates 
with square lattices increase the robustness compared to 
typical hexagonal lattices, such as graphite. Moreover, 
unlike hexagonal lattices, square lattices do not have a 
scan line where the stability becomes extremely weak. 
This has implications for experiments where many dif- 
ferent scan lines are explored, and thus a nanocrystal at 
a stable orientation for one scan line would be forced to 
rotate when the apparatus switches to another scan line, 
as was noted for the specific case of graphite flakes on 
graphite in Ref. [6j. 

Together with the increased moment of inertia for a 
nanocrystal, our results suggest that it should be easier to 
experimentally observe incommensurate sliding for three- 
dimensional nanocrystals with bec or fee structure than 
for hexagonal graphite flakes. 
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